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ARTICLE INFO ABSTRACT 


Keywords: The effect of carrier gas on nanomechanical properties of GaN thin films grown on (0001) sapphire substrates by 

GaN metal-organic chemical vapor deposition (MOCVD) process, was studied by Berkovich nanoindentation. GaN/H2 

Carrier gas and GaN/Nz films were deposited using hydrogen and nitrogen, as carrier gases, respectively. New method was 

rat aa developed to compute the contact stiffness versus the penetration depth using a single load—displacement curve. 

Nanoindentaticn It was found that the hardness and Young’s modulus of GaN/Hz are lower than those of GaN/N2. Moreover, pop- 

Loading rate in events were revealed for GaN/Hg, in contrast for GaN/N3. Fracture was only manifested in the imprint of GaN/ 
Nz due to its low fracture toughness. Besides, it was disclosed the high sensitivity of loading rate only for GaN/ 
Hp. The GaN/Nz2 sample experiences plastic strain relaxation, while GaN/Hg is a highly stressed sample; hence, it 
has less dislocation density compared to GaN/N2 sample. We demonstrated that the used carrier gas influences 
on GaN epitaxial growth process, which in return deeply affects the resulting dislocation density and dislocation 
plasticity mechanisms. These latter, tailor the nanomechanical properties of GaN samples and the indentation- 
induced deformation behavior. 


1. Introduction 


III-V nitride semiconductors, such as GaN, AlGaN and InGaN have 
attracted considerable attention because of their wide band-gap and 
wavelength in the range of red to ultraviolet, for the fabrication of op- 
toelectronic devices, optical detectors, semiconductor lasers, and other 
broad application spectrum [1,2]. 

In previous decades, intensive researches have been focused on GaN 
thin films in industry and academia laboratories owing to its high 
physical properties and its sustainability with low environmental impact 
[3-5]. Besides, the mechanical characteristics of GaN materials are also 
key factors for their successful application, since for instance, the con- 
tact loading during fabrication or packaging can seriously damage the 
devices’ performance. Likewise, GaN based devices experience large 
thermo-mechanical strains, that affect their properties in service, be- 
sides to growth strains that are inherent during GaN crystal growth and 
heat treatment. The knowledge of GaN mechanical properties and un- 
derstanding the elementary mechanisms responsible of 
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nanoindentation-induced plastic deformation in semiconductors is of 
paramount importance for the fabrication, design and service of GaN- 
based devices [6]. 

In last decades, nanoindentation tests have been widely used as a 
powerful tool for probing the mechanical properties and behavior of 
materials within the atomic scale ranges [7,8]. Oliver and Pharr pro- 
posed an analysis method to compute the hardness and Young’s modulus 
from the load-displacement curve obtained from nanoindentation tests 
using Berkovich indenter [7,9]. During nanoindentation, the 
load—-displacement curve is marked by initial branch attributed to an 
elastic behavior with fully reversible deformation [10]. By continuously 
increasing the load, the onset of the plastic deformation (irreversible 
deformation) appears in general by the observed pop-in event or 
excursion in depth. The first pop-in events were observed on electro- 
polished surfaces of gold, copper and aluminum by Gane and Bowden 
{11]. This phenomenon was also commonly observed in several inden- 
ted semiconductor materials [12-14]. Nevertheless, the threshold of 
plastic deformation is not usually revealed by the appearance of such 
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phenomenon and the elastic-plastic transition occurs unperceived 
without any discontinuities on the load—displacement curve [15]. The 
origin of this difference between the presence of pop-in events and their 
absence is not clearly understood at present. 

Previous researches have reported on the determination of me- 
chanical characteristics of GaN thin films. Tsai et al. [16] reported on 
Berkovich nanoindentation-induced deformation of GaN thin films 
grown using MOCVD process on a-plane sapphire substrate and 
demonstrated that pop-in event is related to the formation and propa- 
gation of dislocations, rather than by indentation-induced phase trans- 
formation. Later, Fujikane et al. [17] performed strain rate controlled 
nanoindentation test on bulk GaN crystals with large range of indenta- 
tion load velocity. They investigated the effect of the load rates on the 
incipient plasticity of GaN crystals. 

Even though, physical properties of semiconductors are affected by 
dislocations that are generated in lattice-mismatch during crystal 
growth process [18]; however, scarce works reported on thorough 
interpretation of various mechanical behaviors observed in semi- 
conductors during nanoindentation-induced plastic deformation. 

In this study, Berkovich nanoindentation tests, were performed on 
two GaN thin films grown on (0001) sapphire substrate by MOCVD 
under similar growth conditions, but with dissimilar carrier gases; H2 for 
the first sample and Nz for the second one. The carrier gas type for the 
GaN growth has attracted many interests [18,19], Amano et al. [20] 
initially adopted pure Hz as the carrier gas. Later, Cho et al. [21] 
adjusted the carrier gas from Hz to pure N2 for GaN growth. The effect of 
carrier gas was studied on the structural, morphological, optical and 
electrical GaN epilayer characteristics [2,21,22]. 

Although of these interests in the carrier gas effects on physical 
properties of GaN films, to the best of our knowledge, the specific 
dislocation mechanisms that tailor the nanomechanical behavior and 
properties associated to the impact of carrier gas are not yet clearly 
addressed and well understood. In this framework, the present research 
provides a systematic investigation of the carrier gas Hz and Ng effect on 
nanomechanical properties and behavior of wurtzite GaN thin films 
grown by MOCVD. We attempt to shed some lights on the relations 
between the physical and mechanical features of GaN epilayers to un- 
derstand their perceived dissimilar nanomechanical behaviors. This 
endeavor is a fundamental challenge, because it might contribute to 
deeply understand the nanoindentation-induced deformation mecha- 
nisms and then offers to fabricated GaN films the target physical and 
nanomechanical properties, leading to achieve high-performance opto- 
electronic devices based on nitrides. 


2. Experimental details 


GaN samples were prepared either with hydrogen (Hz) or with ni- 
trogen (Nz) as carrier gas during growth process at 1100 °C. These GaN 
films were grown on (0001) c-plane sapphire substrates by metalorganic 
chemical vapor deposition (MOCVD) technique. Trimethylgallium 
(TMG) and ammonia (NHs) were only used as precursors of gallium and 
nitrogen, respectively. During the growth process, and for each sample, 
a buffer layer of GaN (30 nm) have been deposited using a different 
carrier gas (either Hz or Nz) under sapphire substrate. The thickness 
attained for both GaN samples is 2 um. To seek for simplicity of the 
presentation, more details regarding the growth conditions are reported 
elsewhere [2,22]. 

In order to analyze the residual stresses within both GaN samples, 
micro-Raman spectroscopy (RS) and photoluminescence (PL) were 
performed. The used PL excitation source is the 325 nm (3.81 eV) 
wavelength line of an He-Cd laser of 10 mW power. Structural properties 
were investigated and dislocation density was revealed using HRXRD 
technique. Hall effect measurements (the standard Van der Pauw 
configuration) were used to obtain the electric properties and to attempt 
to correlate them with the presence of defects that influence the nano- 
mechanical properties of GaN samples. The surface morphology of both 
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GaN samples was assessed using atomic force microscopy (AFM) 
technique. 

Nanoindentation tests were performed using nanoindentation in- 
strument (NanoTest NT1, NanoMaterials, Ltd) equipped with a diamond 
pyramid-shaped Berkovich-type indenter. The three-sided pyramidal 
Berkovich shape has a half angle of 65.3°. The nanoindentation system 
was fully calibrated using fused silica as a standard sample prior to the 
measurements. 

The tip-end curvature radius of Berkovich indenter was also deter- 
mined using Hertzian analysis and found equal to 400 nm. Nano- 
indentation tests were performed at room temperature and atmospheric 
pressure, in the load rate-control mode up to the maximum loads P,,,x of 
10 mN and 15 mN. At least sixteen independent tests were carried out on 
each sample to assess the reproducibility of the tests. At the maximum 
load, the indenter was maintained at a constant load for a dwell period 
of 30 s for the thermal drift reason. Each indentation was separated by 
30 ym to avoid possible interferences between neighboring indents. It is 
worth to mention that all indentation depths are less than 10% of the 
film thickness, where the sapphire substrate effect on the characterized 
properties of GaN films is negligible. In addition, scanning electron 
microscopy (SEM) was used to observe and analysis the indentation 
imprints for both GaN samples. Using Vickers microhardness test, the 
length of generated cracks was measured to evaluate the fracture 
toughness Kjc. 

Ina final step, nanoindentation tests were carried out at four loading 
rates, 0.3 mN/s, 0.5 mN/s, 1 mN/s and 3 mN/s, to study the mechanical 
behavior of GaN samples undergoing an increase of loading rates. 


3. Nanoindentation test analysis method 


Hardness and Young’s modulus are the most frequently mechanical 
properties extracted from nanoindentation testing. Oliver-Pharr’s 
method is the most widely used one to extract these properties from load 
— displacement curves [9]. 

Hardness is defined as the ratio of the maximum load to the projected 
area of the indentation, denoted by A;: 

Prnax 


a a 


c 


where, Pmax is the maximum load and A, is the projected area. 
For a perfect Berkovich indenter, the projected area (contact area) is 
given by: 


Ac = 24.5 he (2) 


where h, is the contact depth and can be calculated as follows: 


Pinax 
he = Nnax — € << (3) 
where, hmax is maximum depth, S is the contact stiffness and ¢ is a 
constant depending on the geometry of indenter. For Berkovich indenter 
tip e = 0.75. 
The reduced elastic modulus of the film is dependent on the contact 
stiffness S and on the projected area, then it is defined as follows: 


1 1 
E, = ap re (4) 


where, S = dP/ 4, is the slope of unloading curve at the maximum load 


point and / is the shape constant equal to 1.034 for Berkovich indenter. 
The relationship between the elastic modulus of the tested film, the 
indenter and the reduced modulus is expressed by: 


1 1-v 1-¥# 


E,  E E; (5) 


where (E, v) and (E;, v;) are the Young’s modulus and Poisson’s ratio 
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of the film and the indenter, respectively. For Berkovich diamond 
indenter, the Young’s modulus and Poisson’s ratio are given respectively 
byE; = 1141 GPa andy; = 0.07 [7]. 

Oliver and Pharr [7] found that the unloading curve is usually not 
linear, and it can be fitted using a power-law relation: 


P(h) = a(h— hy)" (6) 


where, P is the contact force, his the indenter depth, h; is the residual 
depth after indenter withdrawn, and a and m are power-law fitting 
constants. 

Fig. 1 displays a typical load—-displacement curve of nanoindentation 
test, where the contact stiffness S is shown. In certain conditions, it is 
interesting to assess not only the Young’s modulus and hardness at a 
given depth, but also to evaluate the mechanical properties versus 
indenter depth. Indeed, the obtained results are rich of information and 
display the local plastic deformation behavior. Thus, in order to measure 
the nanomechanical properties as a function of the indenter penetration, 
a set of nanoindentation tests should be carried out at several indenta- 
tion stepwise depths. As an alternative to this onerous method, the 
continuous stiffness measurement (CSM) emerges as technique to 
continuously measure the dynamic contact stiffness during the nano- 
indentation test. However, this technique is not always available and 
needs somehow fine adjustment [23]. Instead, in this work, we devel- 
oped a new approach to compute the contact stiffness continuously as a 
function of the penetration depth using only a single load-displacement 
indentation curve performed for a given maximal load (Pmqx). Therefore, 
the mechanical properties are derived as a function of the depth. 

In this new method, starting from a single measured load-unload 
displacement curve obtained at a maximal load Pax. The measured 
load—-displacement curve is subdivided into k-intervals of Ah length, as 
similar as a multiple load-displacement curves obtained at increased 
load P; (see, Fig. 1). This method is based on the assumption of a purely 
elastic unloading path, which is properly justified for some kind of 
materials, where their mechanical behavior is time-independent. In this 
case, the unloading curve is fitted by the power-type law given in Eq. (6), 
as can be suggested by Oliver and Pharr [7]. In order to corroborate the 
elastic unloading path assumption, multi-cyclic nanoindentation was 


P (mN) 


Si 


*~ 


”N 
xz 


h (nm) 
hy bj by Byax 
Fig. 1. Schematic illustration of typical load—depth curve obtained in nano- 


indentation and the definition of the contact stiffness at different loading 
increment based on the assumption of purely elastic unloading. 
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carried out at a maximal load of 15 mN on both GaN samples. One can 
see that for five achieved cycles, the loading path of the next cycle 
matches with the unloading path of the previous cycle, as shown in 
Fig. 2. These results provide a very strong evidence that this kind of 
materials obey a purely elastic unloading path. 

The derived unloading curve corresponding to the load P; is the 
branch of the measured unloading curve translated by a depth increment 
Ah to intercept the measured loading at the P; load value. For each load 
value Pj, the contact stiffness S; is the slope obtained at the maximal load 
P;; where (i = 1..k) corresponds to the number of the increment loads. As 
the a,m and hy parameters are initially determined by fitting the 
measured unloading curve, the contact stiffness S; is calculated at the 
point (Pj, hj + Ah) on the measured unloading curve using the following 
relation S$; = am((hj + Ah) — hy)". This contact stiffness S; is assigned 
to the point (P;, hj) on the unloading curve (i.e., the branch of the 
measured unloading curve), where P; and h; are the maximal load and 
depth at the i” point of the measured loading curve. This calculation 
procedure is applied to determine the different contact stiffness at 
different depth points along the loading curve. This is essential to 
calculate the contact depth h, (see, Eq. (3)), then to obtain the contact 
area A,, which enables to determine the harness H and the reduced 
Young’s modulus E,. 

Commonly, indenters used in nanoindentation are imperfect by na- 
ture, that means the indenter tip feature is not perfectly sharp and has an 
imperfection assimilated to a round end defined by a tip radius. Hence, 
the Eq. (2) expresses the projected contact area along with contact depth 
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Fig. 2. Plots of experimental multi-cyclic nanoindentation load- displacement 
curves obtained with Berkovich indenter at a peak load of 15 mN: a) GaN/H; 
b) GaN/Np. The loading path of the next cycle is overlying the unloading path of 
the previous cycle, illustrating a purely elastic unloading behavior. 
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for the Berkovich indenter is only a theoretical assumption, which un- 
derestimates the effective contact area. Actually, to take account of the 
indenter tip imperfection in the determination of the projected contact 
area, Oliver and Pharr have proposed an empirical formula, given as 
follows [7]: 


Ac = f (he) = 24.502 + Cyhe + Coht!? +... + Cght! 8 (7) 
where, C; through Cg are the fitting constants obtained during the 
calibration process using fused silica with perfectly known mechanical 
properties. 


4. Results and discussion 
4.1. Physical properties of GaN samples 


Commonly, the buffer layer deposited on sapphire substrate at low 
temperature (~600 °C) undergoes annealing stage, where the temper- 
ature increases to prepare the main GaN growth process at high tem- 
perature (~1100 °C). During buffer layer annealing, the surface 
morphology of the buffer layer is affected by carrier gas, among other 
factors [18]. For instance, using Hz carrier gas, the GaN buffer layer was 
transformed into low density isolated nuclei, which were etched by Hz at 


0.200 pm/div 
Z 150.000 nm/div 


1 X 1.000 pm/div 
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hot temperature (see, Fig. 3a). The etched GaN buffer layer was partially 
scattered to form large size GaN islands [24]. These islands have grown 
larger and spread out to coalesce with neighboring larger islands during 
the start-up of the hot temperature GaN growth. This growth regime is 
commonly designated by a three-dimensional growth mode (3D), which 
was observed for GaN grown with H2 carrier gas [22,25]. In this case, 
relatively long duration was necessary for the islands to expanse later- 
ally and to swap the 3D growth mode by a complete 2D lateral growth 
regime, which subsequently has led to a smoother surface, as shown in 
Fig. 3c. 

In contrast, using N2 carrier gas, the produced density of nuclei in 
GaN buffer layer was higher with reduced nuclei size during the 
annealing process, as shown in Fig. 3b. The GaN nuclei coalescence was 
closely immediate owing to this narrow neighborhood of small islands 
that were narrowly close. This 2D growth mode was achieved layer by 
layer along the total growth process of GaN/Nz2 epilayer at hot tem- 
perature (see, Fig. 3d). 

Commonly, the analysis of AFM data gives to quantitative informa- 
tion on the surface roughness. Root-mean square roughness Rpms is 
defined as standard deviation of the surface height profile from the 
average height. Rams was calculated using the following equation: 


X 0.200 pm/aiv 
Z 50.000 nm/aiv 


X 1.000 pm/div 
2 100.000 nm/div 


Fig. 3. AFM images of GaN buffer layer after annealing at 1100 °C with: a) He carrier gas; b) Nz carrier gas. AFM images of GaN epilayers: c) GaN/H2 film; d) GaN/ 


Np film. 
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(8) 


where, N is the number of data points, z; and Z are the height of i data 
point and the average height of the image, respectively. 

For GaN/H2 sample the Rpmys roughness value is equal to 0.5 nm and 
for GaN/N2, Rams roughness was found equal to about 5 nm. 

Indeed, the coalescence time is greatly reduced compared to the 
growth of GaN/H2 sample [2]. This morphological difference between 
the buffer layers grown with Hz and with N2 was attributed to the dif- 
ference in surface energy. 

In order to assess the structural quality of both GaN films, it is 
necessary to compare the intensity and the full width at half maximum 
(FWHM) of their X-ray rocking curves. Fig. 4 a and b display the sym- 
metrical (002) and asymmetrical (105) X-ray rocking curves of GaN/H2 
and GaN/Ng. These results indicate the high crystal quality of GaN/H2 
sample owing to its higher intensity and lower FWHM compared to 
GaN/N2 sample. In fact, the measured rocking curve FWHM values of 
(002) and (105) are 259 arcsec and 262 arcsec for the sample grown in 
H2 carrier gas atmosphere and 830 arcsec and 896 arcsec for Nz atmo- 
sphere, respectively. This is often related to defects, mainly the density 
and distribution of threading dislocations (TDs) developed during the 
coalescence of the GaN growth process and mostly depending on the 
initial deposition conditions. The measured rocking curve FWHM for 
(hkD reflection is given as follows [26]: 


B? (hkl) = p2.(hkl) + f2 (khl) (9) 


where, /,(hkl) is the rocking curve broadening caused by angular rota- 
tion at dislocations and 3, (hkl) is the rocking curve broadening caused 
by the strains surrounding dislocations. Eqs. (10) and (11) relate the 
measured FWHM rocking curves to the dislocation density. 


f(hkl) = 2nIn2b’D = Ky (10) 


B2(hkl) = K.tan°0 (11) 


where, b is the length of the Burger vector and D is the dislocation 
density. Fig. 4c displays a straight line which define a linear relationship 
between the /?(hkl) and the tan?0 for five rocking curves at different 
value of angle 6, where the intercept is K, and the slope is K,. The edge 
and screw dislocation densities were determined using Eqs. (10) and 
(11), respectively. The threading dislocation density is the sum of edge 
and screw type dislocation, for each one corresponds a Burger vector 
(Dedge = 0.3186 nm and Dgcrew = 0.5185 nm). In this context, it was 
found that the threading dislocation density is 5 x 10° cm~? for GaN/N2 
and is equal to 2 x 10° cm~? for GaN/H2. These defect amounts are in 
accordance with the electron carrier mobility that are equal to 250 cm?/ 
V.s for GaN/H2 sample and equal to 102 cm?/V.s for GaN/N2 sample, 
which are related to less defects in the GaN/H2 epilayer compared to 
GaN/N2. 

GaN large isolated islands that were formed in the initial growth 
using Hz carrier gas were found to be truncated shapes, on which the 
lateral growth was achieved during the coalescence of these large 
islands [18,25]. The threading dislocations likely blocked at the trun- 
cated facets and quench; consequently, reducing the propagation and 
the density of threading dislocations. It is worth to note that during the 
initial growth and along the annealing step, the formation of large 
islands was disseminated and etched by the Hz carrier gas, as mentioned 
above; so, the distribution of threading dislocations throughout sample’s 
surface and thickness may undergo an inhomogeneous redistribution. 

In contrast at the initial growth step, the GaN/N2 epilayer was 
featured by the formation of the small islands, in which the threading 
dislocations emerge at the layer surface, as a consequence of the strain 
relaxation effect and leading to high threading dislocation density in the 
sample. This statement was confirmed by the obtained dislocation 
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Fig. 4. a) Symmetric (002) rocking curve reflection; b) asymmetric (105) 
rocking curve reflection; c) square of rocking curve line width / (hkl) versus 
tan?(0) for GaN/H2 and GaN/N2 samples. 


density value [18]. 

Therefore, it is clear that the used carrier gas dictates the incipient 
density and distribution of dislocations and defects resulting from lattice 
mismatch of GaN, deposited on sapphire substrate [27,28] at the early 
stage of the growth (MOCVD) process. 

An earlier report points out that the coalescence phase strongly in- 
fluences the strain relaxation in the layer by inducing high threading 
dislocation density as well as the point defects [29]. 

In order to assess the stress states in GaN films grown using different 
carrier gas, actually micro-Raman spectroscopy is the suitable technique 
enabling this valuation. RS allowed to determine the crystalline quality 
and strain relaxation in the layer, which is directly related to the 
threading dislocation density appraised from FWHM rocking curves. The 
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intense Ez phonon frequency mode [30,31] was used as reference peak 
to determine and compare the compressive biaxial stresses in both GaN 
samples. Fig. 5a shows Raman shift of the E: phonon mode for both 
GaN/Hg2 and GaN/Nz2 epilayers. It is easy to observe that Eo) phonon 
mode frequency for GaN elaborated with Hz carrier gas is located at 
570.04 cm~ and the Ex") phonon mode frequency for GaN grown with 
Nz atmosphere is positioned at 568.13 cm~!. The dashed line refers to 
the strain-free frequency of the Ez phonon mode for GaN bulk layer, 
which is situated at 566.34 cm~! [21]. 

Based on these results, it is important to highlight that the two layers 
are subjected to compressive residual stresses, because they are on the 
right side of the bulk GaN position. 

Moreover, these results show that the GaN/N> film is less compressed 
than that of the GaN/Hg film. Indeed, the strain relaxation of the GaN/ 
N2 film is related to the increase of the threading dislocations embedded 
in the sample, as mentioned in previous reports [32]. In other words, the 
GaN/Hz film experienced less strain relaxation resulting in less thread- 
ing dislocation density, as confirmed by the measured density disloca- 
tion for both GaN samples. The thermal expansion coefficients and 
lattice mismatch at the film and substrate interface are generally the 
principal origins in the difference of the strain relaxation between the 
two GaN samples. Moreover, it appears that the strain relaxation of the 
epilayers is inhomogeneous across the thickness of samples. Indeed, the 
compressive residual stresses are smaller close to the substrate than to 
near the GaN film surface, as reported by Cho et al. [21]. 

Besides, the structural properties of GaN thin films related to the 
residual stress was investigated using PL at temperature of 5 K by 
identifying the excitonic line position. Fig. 5b shows for GaN/H2 higher 
PL peak intensity, with the excitonic line position of 3.483 eV, compared 
to that obtained for GaN/N2 (3.479 eV). The dashed line refers to the 
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Fig. 5. a) Raman spectra at room temperature; b) PL spectra measured at 5 K 
obtained for both GaN/H2 and GaN/Nz epilayers. 
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strain-free excitonic line position for GaN bulk layer, which is placed at 
3.469 eV [33]. The difference in the peak PL position is related again to 
the dissimilar compressive residual stress states for both GaN samples. 
As known, for GaN grown on sapphire, the excitonic lines shift toward 
higher energies as the residual stress in the layer increase [34]. More- 
over, the FWHM of the excitonic peak for GaN/Hg2 epilayer (FWHM (Ha) 
= 4 meV) is smaller than for the other sample (FWHM (Nz) = 10 meV). 
Once again, these results confirm the superiority of the crystalline 
quality of GaN/H2 epilayer compared to GaN/N2 sample. Consequently, 
as above stated, it is clearly shown that both Raman and PL results are in 
well accordance; hence, the crystalline quality and the stress relaxation 
of both GaN samples were efficiently established. Moreover, these re- 
sults are in good correlation with the measured threading dislocation 
density obtained by rocking curve reflections. These findings are in good 
agreement with literature reports [35]. 

Currently, the key question is what are the consequences of the used 
carrier gas on the nanomechanical behavior and properties of GaN 
epilayers? In the following, we show the nanomechanical results and 
based on the analysis of the indentation induced-deformation curves, we 
attempt to interpret the nanomechanical features and behavior in rela- 
tion with the physical properties of both GaN samples. 


4.2. Nanomechanical properties of GaN samples 


Fig. 6a and b show the loading-displacement (P-h) curves for GaN 
films, grown using Nz (GaN/N2) and Hz (GaN/Hg) as carrier gases, 
respectively. Indentation experiments were carried out at two maximal 
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Fig. 6. Load-displacement curves of 16 nanoindentation tests obtained at 
maximal load of 10 mN and 15 mN for: a) GaN/Hz film (red curves); b) GaN/N2 
film (blue curves). 
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loads, 10 mN and 15 mN. The indenter was withdrawn using the same 
loading rate of 0.3 mN/s as in the loading step. It is worth to note, that 
for both samples, the 16 nanoindentation curves are overlapped 
demonstrating the reproducibility of the experiments. 

By analyzing Fig. 6a and b, one can observe that the loading curves 
for GaN samples are different. Fig. 6a displays the occurrence of the pop- 
in event for GaN film grown using Hz carrier gas. This first pop-in event, 
corresponds to a sudden burst during the loading curve, which occurs at 
low load, designed by critical load ofP,, = 1.48 mN (pop-in plateau) fora 
penetration depth of about 40 nm. The pop-in length of the horizontal 
plateau is about 7.2 nm. The dissipated energy for the creation of this 
pop-in event was evaluated at 8.3 x 10~!J. It may be noted that pre- 
vious results reported by Jian et al. [36] are in agreement with the 
current findings. Moreover, at peak loads (i.e., 10 mN or 15 mN) and 
during the dwell period (30 s), the indenter continued to penetrate into 
both samples. The aim of holding the peak load for a sufficient period of 
time was to avoid any effect of the thermal drift on the computed 
hardness and Young’s modulus and to totally guarantee the material’s 
elastic recovery. Additionally, one can observe that due to its local 
resistance to plastic deformation, the penetration depth for GaN/Nz2 film 
was lower than that observed for the GaN/Hp. 

Based on preliminary analysis of these results, it is more likely to 
mention that the usage of Hz gas for the elaboration of GaN films con- 
tributes to a certain ductility and easiness to penetrate into the sample. 
In contrast, the usage of Nz carrier gas to grow GaN films, enhances the 
indenter penetration resistance into the sample. These observations 
were corroborated by the calculated plastic work W, for both GaN 
samples at a maximal load of 10 mN. The calculated plastic works are 
(368 + 10)x107!7J and (265 + 10) x10~!2J for GaN/H2 and GaN/No, 
respectively. This difference in the plastic work reveals the enhanced 
ductility of GaN/H2 film compared to GaN/No. Thus, it is of paramount 
importance to understand the foundations behind such distinct behavior 
of GaN samples. Next, based on thorough analysis of qualitative and 
quantitative results, further interpretations of the observed distinguish 
nanomechanical behavior of the elaborated GaN samples would be 
possible. 

The usage of the developed method permitted to compute the 
Young’s modulus E and the hardness H as function of penetration depth. 
Our chief aim by measuring the depth-dependent hardness of the GaN 
samples was to examine the presence or not of this dependence associ- 
ated to the indentation size effect (ISE). 

Fig. 7a and b show for both GaN samples, the evolution of these 
nanomechanical properties obtained from the average curve represent- 
ing the 16 experimental indentation tests carried out at maximal load of 
10 mN. As can be seen, depth-dependent Young’s modulus and hardness 
are larger for GaN/N2 than that for the GaN/H2. Fig. 7a shows Young’s 
modulus for GaN/N2, which was observed continuously decreasing 
along with the penetration depth. In contrast, for the GaN/Hz sample, 
the evolution is practically constant around an average value of 310 GPa 
and almost independent of the magnitude of the local plastic deforma- 
tion in the range of penetration depth of 50 nm to 150 nm. The peak 
depicted in these curves was associated to the pop-in event observed for 
this sample. Fig. 7b shows the evolution of the hardness as a function of 
the penetration depth for both GaN samples. One can observe the con- 
stancy of the hardness along with the penetration depth after the initi- 
ation of the pop-in event (i.e. beyond almost 40 nm of depth), for the 
GaN/H2. However, for the GaN/Ng, it can be seen the gradual increase of 
the hardness as a function of the penetration depth. This latter increases 
started after a depth of about 30 nm, which indicates the onset of the 
plastic deformation. 

The incipient plasticity is attributed to a maximal shear stress at 
which plastic deformation starts locally beneath the indenter tip. Using 
the Johnson’s model [37], the maximal shear stress was calculated and 
given by the following equation: 
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Fig. 7. a) Young’s modulus; b) hardness as a function of penetration depth of 
GaN/H2 and GaN/Ng2. Black square and disk solid symbols refer to Young’s 
modulus and hardness, measured based on single load-displacement curves 
obtained at various maximum loads. 
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where Po, , E and R are the critical indentation load, the Young’s 
modulus and the radius of the indenter tip, respectively. This critical 
load corresponds to the maximal load at which the local plastic defor- 
mation is initiated underneath the indenter tip. The calculated maximal 
shear stresses induced-plasticity are 17.3 + 0.6 GPa and 16.2 + 0.6 GPa 
for GaN/H2 and GaN/Ng, respectively. 

At the maximal load of 10 mN, the average values of Young’s 
modulus E, the Coulomb’s or shear modulus G and the hardness H are 
313 + 6 GPa, 125 + 2 GPa and 13.7 + 0.6 GPa for GaN/H2 and 356 + 4 
GPa, 142 + 3 GPa and 19.5 + 0.4 GPa for GaN/Ng, respectively. 

As Taylor hardening model represents a straightforward relationship 
between the dislocation density and the indentation hardness, since the 
hardness is proportional to the square root of the dislocation density 
[38]. Consequently, the measured hardness was found in accordance 
with Taylor model prediction because the dislocation density found for 
the GaN/Nz2 epilayer is higher than for the GaN/Hg film. 

In addition, Fig. 7 shows the plotted values of Young’s modulus and 
hardness (Black square and disk solid symbols) obtained from the multi- 
cycle indentation tests. One can see the good agreement between the 
provided results using the new method, developed to calculated the 
stiffness and then to determine the nanomechanical properties as a 
function of the indenter depth, and those obtained from single 
load-displacement curves. This represents a quantitative validation of 
the proposed method to calculate continuously the evolution of 
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nanomechanical properties along with penetration depth using a single 
nanoindentation test. 

Table 1 presents nanomechanical properties of GaN thin films re- 
ported in previous studies of literature and the results obtained in this 
study. It clearly depicts discrepancies between the reported nano- 
mechanical properties for GaN films. This scatter of these nano- 
mechanical properties obtained by various nanoindentation tests is 
mainly attributed to the specific tip — surface contact configuration and 
to the different stress distribution inherent into specimens. This is 
somewhat owed to elaboration conditions, to defect interactions with 
indenter tip inducing local plastic deformation and to the used indenter 
tip-type [39]. 


4.3. Dislocation mechanism scenarios correlated to the epitaxial growth 
processes 


Attaining this stage in the investigation, we attempt to thoroughly 
analysis the correlation between the experimental nanomechanical 
properties and the type of carrier gas used in the GaN growth process. As 
abovementioned, GaN/H2 sample showed a first pop-in event, while for 
the GaN/N2 sample, no pop-in events were observed along the loading 
curve. In fact, several possible interpretations for the manifestation of 
pop-ins were reported in literature. The occurrence of pop-ins was 
sometimes related to phase transformation [44], initiation and propa- 
gation of micro-crack [45] and initiation of plastic deformation, marked 
by the transition between elastic to plastic of the strained locally volume 
underneath the indenter [46,47]. In the present paper, the latter 
explanation is very prospected, as it attributes pop-in events to a 
mechanism of homogeneous dislocation nucleation of loops. Micro- 
scopical investigations have revealed the homogeneous nucleation of 
dislocations process during nanoindentation [48-50]. This dislocation 
mechanism was also strongly supported using molecular dynamic sim- 
ulations (MDS) defect free-volume subjected to critical shear stress, 
which activates slip systems and then induces the nucleation and 
emission of dislocation loops [51]. In other words, homogenous dislo- 
cation nucleation mechanism was the only possible scenario if no mobile 
dislocations already existed in the sample close to the indenter. This is 
likely encountered in some materials featured with lower pre-existing 
mobile dislocation density that were not emerging at the sample sur- 
face. Previous works reported that many factors influence on the me- 
chanical behavior and particularly on the initiation of plasticity in thin 
films [44,52,53]. Indeed, the crystallographic texture, the presence of 
surface roughness [54,55], interstitial atoms, density and distribution of 
pre-existing dislocations [6] and the influence of indenter tip geometry 
{39] are among the aspects that influence the appearance or disap- 
pearance of the pop-ins in material indentation. According to literature, 
dissolved interstitial atoms in crystalline materials, featured by shallow 
dislocation density, influence the pop-in event occurrence [50]. More- 
over, the pop-in event positions were also strongly influenced by the 
interaction between the pre-existing threading dislocation density and 
the interstitials, where these latter pinned and locked of pre-existing 
mobile dislocations [49,56]. Unlocking of pre-existing dislocations 
that were pinned by interstitials can produce also the appearance of pop- 


Table 1 
Nanomechanical properties of GaN thin films compiled from literature with 
findings of the present work. 


References E(GPa) H (GPa) Tmax (GPa) 

Jian [36] 314.93 + 40.58 19.34 + 2.13 - 

Tsai [16] 286 + 25 1941 6.3 

Chien [40] 286.12 + 25.34 19.31 + 1.05 - 

Kucheyev [41] 233 13.4 11.5 

Fujikane [42] 323.8 20.0 18.941.3 

Huang [43] 333.6 + 2.7 19.04 + 0.23 - 

This work GaN/N2 356 + 4 19.5 + 0.4 16.2 + 0.6 
GaN/H2 313 +6 13.7 + 0.6 17.3 + 0.6 
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ins [57]. 

In this context, for a better argument foundation of the observed pop- 
ins phenomenon in GaN/H2 sample, we attempt to evaluate the density 
of dislocation loops using the classical dislocation theory. Actually, the 
free energy U required for homogeneous emission of circular dislocation 
loop of radius r and Burgers vector magnitude b, by applying the shear 
stress Tmax to initiate plastic deformation, is given by the following 
equation: 


U = 297 tig — TmaxbAV" (13) 
where 74, is the line energy of the generated dislocation loop. Indeed, 
the Eq. (13) is composed of two terms. The first term expresses the 
necessary energy to create one dislocation loop in free-defect lattice and 
corresponds to the increased lattice energy due to the nucleation of the 
circular loop. The second term represents the necessary energy as the 
external work for the expansion of the dislocation loop. The line energy 
Yais is given by the following equation: 


tas = G = “| ° (we 2) (14) 
where, G, v and ro are the shear modulus, Poisson’s ratio and the 
radius of the dislocation core of GaN/Hz film, respectively. On one hand, 
the free energy U has a maximal value, when the dislocation loop radius 
attain a critical value r, and on the other hand, it is in the similar order of 
the thermal energy kT ~ 10~?J (k is the Boltzmann constant and T is the 
absolute room temperature) which is very small, U(r.) ~ 0. Thus, by 
considering these conditions, the following relations are derived: 


andr, = ee). 
bt max 4 


27, 
fos Ydis 


c 


(15) 


Using these equations, the estimated critical and core radii for the 
dislocation loop arer, ~ 1.8 nm and ry ~ 0.36 nm for the GaN/Hg film, 
respectively. The number of dislocation loops generated at the first pop- 
in is estimated using the ratio between the free energy which is esti- 
mated around of U ~ 7 x 10717 J and the dissipated energy for the 
creation of this pop-in event (evaluated at 8.3 x 107!7J). The obtained 
number of dislocation loops associated to the pop-in event is ~ 10° 
dislocation loops. This number is considered small and it is in accor- 
dance with the homogeneous dislocation nucleation mechanism [58], 
rather than with the collection activation of pre-existing mobile dislo- 
cation scenario (Dislocation cross-slip in GaN single crystal). Using this 
method and assuming that the total dissipation energy in the nano- 
indentation was completely transferred into the formation of dislocation 
loops with critical radius, the number of the formed dislocation loops is 
in the order of 5 x 10° dislocation loops. This value is regarded as an 
overestimation of what can be the upper limit of the number of the 
formed dislocation loops, when the total dissipation energy representing 
the area between the loading and unloading curves was used. 

The GaN/N2 indentation loading -displacement curve shows a 
different nanomechanical behavior of local volume deformed by shear 
stress underneath the indenter tip, compared to that observed for GaN/ 
Hg film. As reported in literature, the pre-existing dislocation density 
usually has a strong effect on the revelation or not of pop-in events 
{13,49,59]. Lorenz et al. [58] stated with a high degree of certainty that, 
the pop-in event will not appear for dislocation density higher than 
10° cm~? for high indenter tip radius. Herein, it is clearly noticed the 
absence of distinct pop-in event (see, Fig. 6b). This is commonly 
attributed to the interaction between the indenter tip and the high pre- 
existing mobile dislocation density. The high dislocation density 
encountered in GaN/N; film is ~ 5 x 10° cm~? and it is directly related 
to the epitaxial growth process using No as carrier gas. 

Consequently, the low or high pre-existing dislocation density 
(determined using X ray rocking curves), was generated and propagated 
since the early growth process of GaN samples. The growth mode (either 
2D or 3D revealed using AFM images) is dictated by the used carrier gas- 


N. Boughrara et al. 


type and therefore affects the residual dislocation density within sam- 
ples. This variation of the pre-existing dislocation density contributes to 
the increase or decrease of the compressive stress, assessed using RS and 
PL techniques, as a consequence of the grade of strain relaxation ach- 
ieved along with the crystal growth. During the nanoindentation test, 
the interaction between the indenter tip and the residual dislocations 
underneath the indenter was manifested by the generation and propa- 
gation of new dislocations according to the related dislocation mecha- 
nism. This latter is expressed on the load nanoindentation curve by the 
presence or not of pop-in events. Indeed, for high residual dislocation 
density, the collective movement and multiplication of mobile disloca- 
tions is the eminent dislocation mechanism-based plastic deformation. 
In contrast, for low residual density dislocation, the induced-plastic 
deformation, underneath the indenter tip, is mostly attributed to the 
homogeneous dislocation nucleation mechanism. 

In order to illustrate the concept of the dissimilar dislocation 
mechanisms encountered in each GaN film during indentation test, 
Fig. 8 shows schematic representations (models) of dislocation mecha- 
nisms that can occur during the nanoindentation of GaN/H2 (Fig. 8a) 
and GaN/N2 (Fig. 8b). The dislocation mechanism is related to the 
incipient dislocation density and to the compressive strain state that 
were originated from the used carrier gas. In this study, we present this 
interpretation as a speculation of the most probable scenarios that might 


a) GaN/H2 
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occur and enable us to understand the distinct nanomechanical behavior 
observed on GaN nanoindentation curves. The speculation of such 
dislocation mechanisms could be visually assessed using transmission 
electron microscopy (TEM) technique as a possibly future work. 


4.4. Fracture toughness assessment 


Nanoindentation test has proven to be a powerful technique to 
evaluate material’s brittleness or resistance to fracture. The fracture 
toughness was evaluated using nanoindentation tests and is defined by 
Lawn-Evans—Marshall (LEM) model as follows [60]. 


E\'? P 
ke =a(5) BR 


where E is the Young’s modulus, H is the hardness, P is the inden- 
tation load, c is the crack length and « is an empirical constant that 
depends on the geometry of the indenter (« = 0.016 for the Vickers 4- 
sided pyramidal indenter). This equation defines the relationship be- 
tween the fracture toughness Kjc and the ratio of hardness to Young’s 
modulus of a material (H/E) [61]. The higher ratio (H/E), the less 
fracture toughness Kjc is. Besides, the H*/E? ratio is defined as suitable 
index to assess the ability of material to resist the plastic deformation of 
the material [62]. 
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Fig. 8. Schematic representation of dislocation mechanisms responsible for the nanomechanical behavior observed in load—displacement curves of: a) GaN/H2; b) 


GaN/No2. 
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Fig. 9 shows the average variation of H/E and H?/E? ratios as a 
function of the penetration depth for GaN/Hz and GaN/Ng films ob- 
tained at the maximum loading value of 10 mN. 

This figure shows that starting from nearly 40 nm of penetration 
depth, H/E and H°/E? ratios obtained for GaN/N2 are higher than those 
obtained for GaN/H2. This means that the fracture toughness Kya, which 
is inversely proportional to the square root of (H/E) ratio of GaN/Nz is 
lower than that of GaN/H2. In other words, the probability of GaN/N2 
film fracture, when critical load of fracture was attained, is higher than 
for GaN/H2 sample. Moreover, the increased penetration depth to pro- 
duce larger plastic strains in GaN/N2 required higher load increment 
compared to GaN/Hg3. In fact, Fig. 9 clearly shows that for GaN/Hp, both 
(H/E) and (H?/E?) ratios are fairly steady around an average value; 
however, for GaN/N2, they are increasing with penetration depth. In 
order to interpret these curve trends, we refer to the evolution of 
Young’s modulus and hardness for both samples, as displayed in Fig. 7. 
For both GaN samples, hardness is around a steady value starting from 
around 50 nm of penetration depth (see, Fig. 7b). 

However, for GaN/Ng the decrease of its Young’s modulus versus the 
penetration depth is higher (larger decreasing of the slope) compared to 
that for GaN/Hz2 (see, Fig. 7a). Based on these findings, GaN/Hg film 
shows a ductile behavior, where it is possible to undergo more plastic 
strains by applying locally sever shear stresses. In contrast, for GaN/N2 
sample, (where the majority of these dislocations emerged at the sample 
surface), the indenter tip interacts with the collective movement and 
multiplication of the activated pre-existing mobile dislocations resulting 
in higher dislocation density underneath the indenter tip. This may lead 
to an increase of the probability of the initiation of internal fracture, 
which can emerge to the sample surface as cracks. Indeed, the decrease 


a) 0.06 

0.04 
= 
= 

0.02 

ie i 1 1 
0 50 100 150 
Penetration depth (nm) 

b) 

0.06 + an 

fe 
& 0.04+ 
gy ail 
a 
a 
” 
= 0.02 + 
—e— GaN/N, 
0.00 £ —e— GaN/H, 
I 1 n 1 
0 50 100 150 


Penetration depth (nm) 


Fig. 9. a) H/E; b) H?/E? variation versus penetration depth calculated for a 
maximum load of 10 mN, for both GaN films. 


10 


Applied Surface Science 579 (2022) 152188 


of the Young’s modulus along with the penetration depth for GaN/Nz is 
synonym to a decrease of the material stiffness, which is attributed more 
likely to the formation of nano cracks that grow and coalesce to lastly 
lead to fracture. In order to quantitatively assess this behavior, the 
calculated fracture toughness was found equal to Kjc = 0.85 MPa m\/2 
for GaN/No2, using LEM model given by Eq. (16). This sample showed 
failure and crack propagations for an applied load of 1962 mN using 
Vickers microhardness indenter (see, Fig. 10a). Moreover, the higher 
residual stresses and their inhomogeneous distribution within the thin 
film can significantly disturb the fracture behavior and scattered cracks 
appears around the imprint [63]. Nonetheless, this finding is compara- 
ble to the value obtained in literature and it is in the range between 0.7 
and 1.0 MPa m2 [64,65]. While in alike conditions, the GaN/H2 loaded 
using Vickers microhardness indenter, no apparent of cracks was shown 
(see, Fig. 10b). However, it looks that a small pile up phenomenon is 
observed around the borders of the imprint, which reflects the large 
capacity of plastic deformation that can undergo the GaN/H2 sample 
and reveals its higher ductile behavior. It is supposed that more loading 
force should be applied to likely attain failure and therefore cracks 
appearance. 

To sum up the main results in this subsection, it is interesting to 
emphasis that crack phenomenon is attributed to the interaction and 
propagation of dislocation activities in the area of the indented surface. 
Whereas, the crack propagation is caused by highly strained field vol- 
ume underneath the indenter and crystallographic properties of the 
material [45,66]. Therefore, the carrier gas used for GaN growth process 
plays a crucial role on the resistance of material to fracture underneath 
an indenter tip or simply to its fracture toughness value. 


a) 


cracks 


cracks 


Fig. 10. SEM images of indented regions with a Vickers microhardness 
indenter at load of 1962 mN on the top surfaces of both GaN films: a) GaN/N2 
sample displays the cracks along corner and edges of the imprint; b) GaN/H2 
sample shows the pile-up phenomenon around the imprint boundaries. 
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4.5. Loading rate effect 


Fig. 11 shows the load-displacement curves obtained at various 
loading rates, namely, 0.3, 0.5, 1, and 3 mN/s for both GaN thin films. 
One can distinguish that load-displacement curves for GaN/Hg2 are 
sensitive to the loading rate. Moreover, the width of the pop-in plateau 
and the critical load are increasing with loading rate (see, Fig. 11a). For 
instance, 6 nm, 13 nm, 21 nm and 29 nm are pop-in plateau width, and 
the associated critical load are 1.48 mN, 4.27 mN, 6.13 mN and 8.46 
mN, obtained for the following loading rates: 0.3 mN/s, 0.5 mN/s, 1 
mN/s and 3 mN/s, respectively. As consequence, for GaN/H2 the shear 
stress is loading rate-dependent. In addition, the work of pop-ins is 
increasing; so, the thermal activation of the homogeneous dislocation 
nucleation mechanism [67,68], which is responsible for the formation of 
loops that contribute to the increase of the density of dislocations, is 
getting higher along with loading rate. Therefore, the local volume of 
material underneath the indenter tip is featured by an increase of the 
hardness along with the loading rate rise. GaN/H2 shows a strain 
hardening effect, which could be explained by the rise of loading branch 
slopes for the load-displacement curves with the increase of the loading 
rate. Such effect was not revealed in the investigated bulk GaN samples 
by Fujikane et al. [17]. 

In contrast, for GaN/N2, the measured load—displacement curves at 
different loading rates are fairly overlapped, even for high loading rates, 
as higher as 3mN/s, (see, Fig. 11b). This assumes that the load associated 
to the onset of plastic deformation is almost the same and then loading 
rate-independent, therefore the maximal shear stress for the transition 
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Fig. 11. Representative load—displacement curves of: a) GaN/Hz2 epilayer, 
where the dynamic effect has an influence on the plastic onset and behavior; b) 
GaN/Ng, where no effect of loading rate was observed on the nanoindentation 
load-unload curves. 
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from elastic to plastic deformation is identical. In this case, it should be 
noted that the increase in the loading rate did not induce any increase of 
the sample hardness, i.e., the plastic deformation at high loading rate 
was induced by rapid movement of collective mobile dislocations that 
probably contributed to an increase of the material damage, due to its 
low fracture toughness Kjc, as previously mentioned. 

Table 2 lists for both GaN samples, the calculated maximum shear 
stress and the average of the plastic work associated to the range of the 
considered loading rates. One can observe the gradual rise of the 
maximum shear stress versus the increase of loading rate, for GaN/H2 
sample. However, the average plastic work is decreasing, which is 
synonym of losing its ductility. This means that the reduction of ductility 
along with the loading rate is attributed to an increase of the sample 
hardness, interpreted as work-hardening effect. In contrast, the 
maximum shear stress and the plastic work are loading-rate independent 
for GaN/N2. 

In order to quantitatively assess the Young’s modulus and hardness 
evolution as function of the loading rate, these results are plotted in 
Fig. 12 for both GaN samples. For GaN/H2, Young’s modulus and 
hardness are increasing along with the loading rate. These results reveal 
its plastic hardening behavior. In the microscopic point of view, this is 
explained by the homogeneous dislocation nucleation mechanisms that 
create and emit dislocation loops in a sample with a low pre-existing 
dislocation density, as depicted in Fig. 8a. However, for GaN/Ng, it 
appears that the Young’s modulus was decreased by increasing the 
loading-rate; that means, the local volume underneath the indenter tip 
was losing its stiffness, as shown in Fig. 12a. While, Fig. 12b shows that 
GaN/N2 hardness was not undergoing a hardening behavior by 
increasing loading rate; but in contrast, its hardness was somehow 
saturated. This behavior is likely interpreted as an obstructive effect of 
the high dislocation density of GaN/Ng, as displayed in Fig. 8b. Thus, the 
possibly dislocation mechanism scenarios, in such behavior, were the 
rapid collective movement of mobile dislocations [69] by increasing the 
loading rate. Accordingly, nano-cracks were generated and emerged at 
the surface’s sample, which reduced the sample stiffness as revealed by 
the decline of the Young’s modulus. 


5. Conclusion 


The present work reports on a comparative study of the nano- 
mechanical behavior and properties of GaN samples elaborated within 
Hz and Nz ambient, as carrier gases using MOCVD process. Nano- 
mechanical properties of two GaN thin films were determined by 
instrumented nanoindentation technique. To derive Young’s modulus 
and hardness depth-dependent from only a single load-displacement 
curve, new method was proposed, as alternative to CSM technique. For 
this purpose, pure elastic unloading path was assumed and verified. 
Depending on the epitaxial growth mechanism of GaN (either 3D or 2D 
growth mode) dictated by the used carrier gas, its significant influence 
on the resulting nanomechanical behavior and physical features was 
revealed. The observed distinguish nanomechanical behavior and 
properties is mainly attributed to the distribution of pre-existing dislo- 
cations and to the dislocation density within the sample as a result of 


Table 2 
Maximum shear stress and plastic work calculated at different loading rates for 
both GaN samples. 


Loading GaN with H2 GaN with N2 
te N, 
- es (mN/ Maximum Average of Maximum Average of 
shear stress plastic work shear stress plastic work 
Tmax (GPa) (10-17J) Tmax (GPa) (10-12J) 
0.3 17.3 368 16.2 265 
0.5 20.6 356 16.2 266 
1.0 24.2 345 16.3 264 
3.0 28.4 327 16.5 266 
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Fig. 12. Nanomechanical properties evolution as a function of the loading rate 
for GaN/H2 and GaN/N2 samples: a) Young’s modulus; b) Hardness. 


either the higher or the lower strain relaxation. Based on classical theory 
of dislocations, the occurrence of the pop-in event is attributed to ho- 
mogeneous dislocation nucleation mechanism for the GaN/H2 sample. 
However, for GaN/Nz no such phenomenon was revealed; thus, the 
collective activation of pre-existing mobile dislocation is the most 
probable scenario. The fracture toughness coefficient Kic was deter- 
mined for GaN/N2, where it was revealed cracks at Vickers microhard- 
ness imprint with a load of 1962 mN. For GaN/Hg, a pile-up was 
observed and no cracks was disclosed at sample’s surface. In addition, 
the effect of loading rate on the nanomechanical behavior were shown 
for GaN/Hg2, in contrary to the GaN/N2 sample, where no effect was 
distinguished. Afterwards, these findings emphasis the decisive effect of 
elaboration conditions for GaN epilayers and permit a clear under- 
standing of the close relationships between the resulting microstructure 
and the relevant changes in the nanomechanical behavior caused by 
specific dislocation mechanisms. 
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